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Abstract 

We present a new experimental method to correlate the isotopic composition of intermediate 
mass fragments (IMF) emitted at mid-rapidity in semi-peripheral collisions with the emission 
timescale: IMFs emitted in the early stage of the reaction show larger values of <N/Z> isospin 
asymmetry, stronger angular anisotropies and reduced odd-even staggering effects in neutron to 
proton ratio <N/Z> distributions than those produced in sequential statistical emission. All these 
effects support the concept of isospin "migration" , that is sensitive to the density gradient between 
participant and quasi-spectator nuclear matter, in the so called neck fragmentation mechanism. By 
comparing the data to a Stochastic Mean Field (SMF) simulation we show that this method gives 
valuable constraints on the symmetry energy term of nuclear equation of state at subsaturation 
densities. An indication emerges for a linear density dependence of the symmetry energy. 

PACS numbers: 25.70Mn, 25.70Pq 
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I. INTRODUCTION 



In the Fermi energy domain (15-100 A.MeV) semi-peripheral heavy-ion collisions are 
characterized by the presence of ternary or quaternary reactions leading to the formation of 
a binary system with excited projectile-like (PLF) and target-like (TLF) fragments strongly 
correlated with one or more IMFs and light particles in the exit channel. The emission 
pattern of these reactions has shown that light charged particles (LCP, Z<2) and IMFs 
(Z>3) are not entirely described by the statistical decay of the PLF and TLF. In particu- 
lar, light IMFs (Z<10) are emitted preferentially towards mid-rapidity region (intermediate 
between PLF and TLF rapidities). Their velocity distributions display typical forward- 
backward asymmetry in the invariant cross section d 2 a/v±dv±dv\\ (see, for example, Fig. 
le) in section II), indicating their dynamical origin jl- 3] . The term "neck fragmentation" 
is commonly used for such a type of events, because a neck-like structure is predicted to be 
formed between the two main residues (PLF and TLF), in the early stages of the reaction. 
Transport model simulations jj] have shown that IMFs in the neck region are formed in 
dilute matter in contact with the PLF and TLF residues. However, the probability to pro- 
duce fragments at mid-rapidity depends on different variables, such as impact parameter, 
beam energy and isospin asymmetry [5]. The effective interactions driving transport phe- 
nomena of neutrons and protons through the neck have been related to the slope (density 
gradient) and magnitude (isospin gradient) of the symmetry potential term of the Equation 
of State (EOS). Thus, heavy ion collisions with projectiles and targets with different isospin 
asymmetries have been studied to probe the density dependence of the symmetry term of 
the EOS ("asy-eos") that is an important ingredient for nuclear structure and astrophysical 
phenomena p, Q] ■ Different observables, mainly based upon measurement of the neutron to 
proton (N/Z) ratio of reactions products have shown sensitivity to the density dependence of 
the simmetry energy like isospin diffusion and equilibration SHlOj], neutron to proton ratio 



measurements 



111 ], transverse collective flow 



12|. 



Neutron enrichment of LCPs and IMFs in the mid-rapidity region has been observed ex- 
perimentally (see for example ref. [ijj and references therein). An excess of neutron (isospin 
migration) towards the dilute low-density region is predicted by transport simulations. This 
is determined by the slope of the density dependence of the symmetry energy, providing a 
drift contribution to the isospin transfer, thus producing neutron rich IMFs at mid-rapidity 
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15] . Such neutron enrichment of the neck-like IMFs has never consistently been used 



16 



171 ] . As shown in [18| the isospin 



to pin down the time scale of the reaction dynamics and to constrain the theoretical simu- 
lations vs. the isospin asymmetry. Indeed, a neutron enrichment of the neck region could 
be also affected by phenomena not strictly related to the isospin dynamics, like proximity 
effects in the decay of PLF and TLF [2j, or effects of reduced size of neck structures and 
the persistence of neutron-rich matter at the surface 
drift and isospin diffusion can simultaneously compete to characterize the midrapidity and 
projectile residue emission. 

In ref 19M22I] a new method based on a three body analysis of fragment-fragment relative 
velocities has been introduced to calibrate the timescale of IMFs emission in semi-peripheral 
collisions, thus probing the dynamics and chronology of fragment formation. 

In this paper, the time scale (chronology) of IMFs originating in ternary reactions is 
correlated with their isospin content by inspecting typical dynamical observable, such as the 
degree of alignment and relative velocity correlations. The results are compared in the last 
section to predictions of microscopic transport calculations for the neutron rich system. 



II. EXPERIMENTAL DETAILS AND RESULTS 



The experiment was performed at the INFN-LNS Super-Conductive Cyclotron of Cata- 
nia (Italy), bombarding thin (« 300/ig/cm 2 ) self-supporting 64 Ni and 58 Ni targets with 35 
A.MeV 124 Sn and 112 Sn beams, respectively. Reaction products were detected with the 
forward part of the 4vr multi-detector CHIMERA [22(| that is constituted by 688 Si(f»300 
//ra)-CsI(Tl) telescopes over a total of 1192, arranged in 18 rings and covering the angular 
range between 1° and 30°. Isotopic identification was obtained up to O xyg en by AE — E 



technique. Details about the experimental methods are discussed in [22|, |23| . In our reverse 
kinematics condition the angular coverage represents almost 85% of the cm. solid angle for 
the reaction channels under study. Complete events (at least 70% of the total charge and 
total parallel momentum of the colliding systems) were analyzed. Semi-peripheral collisions 
were selected gating on the total charged products multiplicity M. Here only events with 
M< 6 will be considered. These events correspond to a reduced impact parameter b re d ~ 
b/b max > 0.7, determined using the Cavata approximation 24| as explained in ref. [25]. As 



shown in Fig. 2 of ref [20], favorable conditions of reverse kinematics and capability of the 




FIG. 1. (Color online) For the 124 Sn+ 64 Ni reaction and M< 6 a) correlations between relative 
velocities V rel/^ Viola of the three biggest fragments of the event for IMFs (3<Z<18; b) distribu- 
tion of cos(9pftox) for Z=3 (dotted line) Z=6 (thick line) and Z=18 (thin line); c) as a) with the 
condition cos(9ppox) >0.8; d) as a) with the condition cos(9ppox) <0; e) Invariant cross-section 
for Z=6 IMFs. f) as e) with the condition cos(9ppox) <0; the circumference shows a Coulomb 
ridge centered at PLF source velocity and 2.1 cm/ns radius. 

CHIMERA device to detect fragments in a broad range of kinetic energies (including the 
slow moving target-like residues) greatly facilitate the distinction of PLF,TLF and IMFs in 
ternary events. 

In order to evaluate the timescale of fragment formation we extended the method quan- 
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titatively described in ref. 19|, |20j. The three biggest fragments in each event were 
sorted according to the decreasing value of their parallel velocity Vn along the beam di- 
rection (l=fastest fragment, 2=slowest fragment, 3=intermediate velocity fragment) and 
the fragment-fragment relative velocities, Vr£l(1,3) and V_&el(2,3), were calculated. The 
fragment labels 1,2 and 3 correspond to the PLF, TLF residues and the IMF fragment, re- 
spectively, mentioned on ref. 20j. In fact, ranking the fragments according to their parallel 
velocity gives the simplest way to extend the present study towards the most dissipative 
collisions, associated with higher values of the IMFs multiplicity [261 ] . The relative velocities 
are normalized to the one corresponding to the Coulomb repulsion, as given by the Viola 
systematics, Vyj ; a (i,3) (i=l,2) |27j. The correlation for fragments with charge 3<Z<18 
between the two relative velocities rj=V i?£z/Vyj j a (l,3) and r 2 =V^£;i/Vyj za(2,3) is shown 
in Fig. la): the values ri=l and r 2 =l correspond to sequential decay of IMFs from a 
PLF and TLF respectively; values of r x and r 2 simultaneously larger than unity indicate 
IMFs of dynamical origin (prompt ternary division). Timescale calibration of Fig. la) was 
done in ref. [li], [3] using a three-body collinear Coulomb trajectory calculation, showing 
a well defined chronology: light IMFs are emitted either on a short timescale (within 50 
fm/c) or sequentially (>120 fm/c) after the re-separation of the binary PLF-TLF system. 
This result has been reproduced by different transport reaction simulations like Stochastic 



Mean Field SMF 



28j | and it is 



4| and CoMD-II Constrained Molecular Dynamics model 
in agreement with simulations based on similar three-body Coulomb trajectory calculations 
29|| . It is interesting to note that antisymmetryzed molecular calculations (AMD) [30|] for 
semi-peripheral reactions also predicts the formation of cluster at mid-velocity on a fast 
timescale, within 90 fm/c, due to the inset of shape and density fluctuations. 

In order to characterize the dynamical emission v.s. isospin degree of freedom, the 
quantity cos(9 PRO x) was evaluated. Oppox gives the angle between the separation axis 
ns = {y c pLF-iMF ~ ^tlf) (relative velocity between TLF and PLF-IMF center of mass) 
and the break-up axis rip = (Vplf — Vjmf) (relative velocity between PLF and IMF ori- 



ented from the light to the heavy fragment), i.e., cos 



(9, 



definition is slightly different from the one used in 31 



PROX) 



^^^j. Notice that this 



l n sl|nF| ' 



32] since it requires the crucial de- 



tection of a TLF fragment, as it is in our case. Recently, contemporary TLF-PLF detection 
at Fermi energies has produced important advances in the fields of mass-energy transport 



phenomena 



33( and decay of very heavy nuclei [34j]. Fig. lb) shows the cos(9ppox) dis- 



FIG. 2. (Color online) For the 124 Sn + 64 Ni reaction and M< 6: <N/Z> for charge Z=6 for different 
bins in the ri-r2 plane. The dashed contour, projected in the ri-r2 plane gives a schematical selection 
of the zone with values of ri and r2 simultaneously larger than unity. 

tribution for Z=3, Z=6 and Z=18 IMFs charges respectively. The strong enhancement of 
the distribution for cos(9ppox) >0.8 as seen in Fig. lb) indicates a clear contribution of 
dynamical emission. Notice the strong tendency (for cos(6ppox) ~1) to a backward IMF 
emission in a strict aligned configuration along the TLF-PLF separation axis. The enhance- 
ment for the Z=18 charge is mainly due to the onset of dynamical fission of the projectile 
as shown in [25]. By setting the condition cos(9prox) <0 (forward emission) in Fig. la) we 
obtain the pattern of Fig. Id), where the events populate the region around the axis where 
ri=l, as expected for a sequential decay of the IMF from the PLF source. In contrast, 
in Fig lc) the correlation of Fig. la) is obtained by selecting events with the conditions 
cos(9ppox) >0.8 (backward emission), showing events concentrated near the diagonal, as 
expected for the dynamical emission of fragments. These results are further illustrated by 
comparing the invariant cross-section distributions of Fig. le) and Fig. If) as a function 
of V||-Vj_ velocities with respect to the beam axis. In Fig. le) the invariant cross-section 
for carbon is shown with no selection in cos(6ppox) and in Fig. If) it is shown with the 



condition cos(9prox) <0. In Fig. If), we observe, consistently with the result of sequential 
emission of Fig. Id), a characteristic pattern of a forward sequential emission from the 
PLF source, well shaped around a Coulomb ring. No such evidence is shown in Fig. le), 
indicating that the PLF backward sequential decay, going towards the mid-velocity region, 
is strongly mixed with a non-sequential prompt emission. 

To get more insights to the correlations between isospin, relative velocities and emission 
time-scale of IMFs we have calculated, for each bin (0.5x0.5 width) in the plane rl-r2 of 
Fig. la), the average N/Z isotopic distributions for all charges between Z=3 and Z=8. Fig. 
2 shows, for example, the result for charge Z=6: the largest values of the neutron to proton 
ratios are obtained for events near the diagonal of the Y REL/^vioiai^fi) vs. V rel /"V VioZa(2,3) 
plane, corresponding to prompt emission (lowest time-scale emission) and to the highest 
degree of alignment. The neutron content enrichment of the mid- velocity emission with 
respect to the N/Z of fragments evaporatedjxom PLF or TLF has been previously observed, 



with different interpretations of the results 



|35H37|. Our experimental approach introduces 



two new aspects: i) the <N/Z> distribution is not averaged over the whole emission time 
scale, but it is linked to the reaction dynamics in a consistent way; ii) the <N/Z> of the 
fragments can be correlated with the alignment properties of different emission scenario, so 
enhancing the experimental sensitivity to select those genuine effects due to isospin dynamics 

We show in Figs. 3a) and 3b) the <N/Z> as a function of the IMFs charge Z for 
the reactions 124 Sn + 64 Ni and 112 Sn + 58 Ni, respectively. Our purpose is to measure the 
degree of neutron enrichment at midvelocity and to compare it with the one related with the 
statistical emission from a PLF source for the two systems with different isospin asymmetry. 
The condition cos(9prox) <0 has been used to select fragments statistically emitted in the 
PLF forward hemisphere (solid squares in Fig. 3). Dynamically emitted IMFs (solid circles 
in Fig. 3) are selected by requiring cos(9prox) >0.8 and by imposing a further condition 
in the rj-rg plane that selects events near diagonal of that plane. This latter condition is 
schematically shown by a dashed line in Fig. 2 for charge Z=6. We clearly observe that the 
N/Z ratio for dynamically emitted particles (DE) shows systematically larger values for both 
systems with respect to the one obtained for statistically emitted particles (SE). A second 
interesting observation is the flattening of the even-odd effect in the <N/Z> distribution of 
the neutron rich system with respect to the neutron poor one and this effect is present also 
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FIG. 3. (Color online) Experimental <N/Z> distribution of IMFs as a function of charge Z for 
statistical emitted particles (solid squares) and dynamical emitted particles (solid circles), for the 
reactions a) 124 Sn + 64 Ni and b) 112 Sn + 54 Ni. 



when comparing the DE distribution respect to the SE one for the neutron poor system. 
Odd-even effects have been linked to the last evaporation step involving just one neutron 



or one proton emission 39M41| or, more generally, to the last steps in the decay channel 



421 ] . Both the explanations are closely linked to the ground state binding energy and level 



density effects of final isotopes 



43 



44j that are mainly responsible for odd-even effects in 



evaporation models. Anyway, as it was stated in |44| and evidenced mainly in Fig. 3b), a 
proper selection of the emissing source shows that odd-even effects are influenced by different 
reaction decay mechanisms. 



III. COMPARISON WITH SMF MODEL 



We have compared our data for the neutron rich system to transport theories using the 
Stochastic Mean Field (SMF) mo del flfl, based on Boltzmann-Nordheim-Vlasov (BNV) 



equation, already used in 



20 



22( to describe the basic experimental features of the PLF, 
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FIG. 4. (Color online) Experimental <N/Z> distribution of IMFs as a function of charge Z for 
statistically emitted particles (solid squares) and dynamical emitted particles (solid circles), for 
the reaction 124 Sn + 64 Ni. Blue hatched area: SMF-GEMINI calculation for dynamical emitted 
particles and asy-stiff parametrization; magenta hatched area: asy-soft parametrization. The 
<N/Z> of primary dynamical emitted IMFs as a function of their atomic number Z, obtained 
from SMF calculations, are plotted in the inset. The hatched zone indicates the error bars in the 
calculations. 

TLF, and IMF in ternary reactions. The SMF model implements the nuclear mean field 
dynamics as well as the effect of fluctuations induced by nucleon-nucleon collisions. Two 
different parameterizations of the potential part of the symmetry term of EOS were used. 
The first one linearly increases with the density (asy-stiff) and the second one (asy-soft) 
exhibits a weak variation around the nuclear saturation density po [ijj . The corresponding 
slope parameter L = 3po(de sym /dp) p=po is around 80 MeV for the asy-stiff and 25 MeV for 
the asy-soft choices, respectively. In the current calculation, isovector thermal fluctuations, 
corresponding to the actual value of the symmetry energy at the neck density, have been 



implemented 46[. Calculations have been performed at 6 fm impact parameter and by 



selecting ternary events, as in the experiment. In order to compare the calculations with 
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FIG. 5. (Color online) For the 124 Sn+ 64 Ni reaction, (empty, black) experimental velocity spectra 
V|| (in the cm. system) of dynamically emitted intermediate mass fragments from charge Z=3 to 
charge Z=6 and (hashed, red) calculated velocity spectra (SMF+GEMINI) for the same reaction. 
Calculations have been normalized to data as explained in the text. The arrows indicate the cm. 
beam velocity as reference. Calculated data are filtered for detector acceptance. 



the data, the primary hot fragments 
using the statistical model GEMINI 



produced by SMF pass throught a de-excitation phase 



471 ] . The average excitation energy of the IMFs before 



the GEMINI step is E*/A^ 2.5 A.MeV. 

Fig. H] shows the calculated <N/Z> as a function of dynamical emitted IMFs atomic 
number. As shown in the inset, the two Asy-EoS parameterizations give rather different 
predictions for the <N/Z> of primary fragments: the asy-stiff parametrization produces 
more neutron rich fragments respect to the asy-soft choice. After the GEMINI secondary- 
decay stage the difference in <N/Z> between the two parameterizations persists for Z<7. 
The asy-stiff parametrization (blue hatched area) matches the experimental data fairly well. 

We have checked if the calculation reproduces some basic features of dynamical emitted 
fragments, like, for example, velocity spectra and charge distributions. In Fig. (5^-d) we 
present the center-of-mass velocity spectra of intermediate mass fragments from charge Z=3 
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to Z=6 selected by requiring the condition cos(9prox) >0.8 and by excluding particles that 
lie along the ri=l and r 2 =l lines. The longitudinal velocity distributions, as observed in 
neck fragmentation, are well centered around the midvelocity region at halfway the velocity 
of the TLF and PLF. In the same figures the experimental data are compared with results 
of SMF+GEMINI calculations for neck emitted fragments in ternary events and asy-stiff 
parametrization (hashed histograms). In the calculated distributions only fragments that are 
originating from a primary "neck" fragment are taken into account. The detector geometry, 
thesholds and time resolutions have been applied to the calculation, simulating in detail the 
response and acceptance of the CHIMERA detector (filter). Calculations were normalized to 
data determining the total yield ratio R = Y(exp)/Y(cal) for all charges between Z=3 and 
Z=9 and assuming the same normalization factor R for each Z. Shapes and relative intensities 
of the experimental longitudinal velocity distributions are remarkably well reproduced by 
the SMF+GEMINI simulation: the calculated distributions, as well the experimental ones, 
extend up to the center of mass velocity in the mid-velocity region; it can be also noticed a 
slight tendency for the experimental spectra to be peaked at higher values of V\\ respect to 
the calculated ones at increasing of the IMF atomic number Z. 

Fig. [6] shows the experimental charge distribution for dynamically emitted IMFs com- 
pared with the calculated charge distribution (square) and filtered by the detector geometry 
(circles). The same normalization factor R of Fig. [5] has been used between experimental 
and calculated data. For charge Z=4, the unbound 8 Be isotope is not included in the data; 



thus it has been filtered also in the calculations. The typical exponential behaviour J4] of 
dynamically emitted fragments charge distribution is fairly well reproduced. 

Fig. EH) (solid circles) shows the correlation between <N/Z> and cos(8ppox) for the 
reaction 124 Sn + 64 Ni. In this case, we have restricted the <N/Z> analysis to IMF charges 
between Z=5 and Z=8, excluding lithium and beryllium, in order to increase the sensitivity 
of the analysis by avoiding <N/Z> strongly reflecting the value of stability line. We observe 
an increase of the <N/Z> for cos(8ppox) values larger than 0.9, corresponding to the highest 
degree of alignment. For comparison Fig. [7b) shows the same correlation obtained when only 
IMFs that lies along ri=l line in the V REi/^vioia plot are taken into account (statistical 
emission). Notice the flat distribution, as expected for a statistical isotropic decay; the 
average <N/Z> value is equal to 1.15. 

In Fig. UJi) the SMF calculations of the <N/Z> values are plotted as a function of 
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FIG. 6. (Color online) For the 124 Sn+ Ni reaction, experimental charge distribution of dynam- 
ically emitted IMFs. The squares correspond to the SMF+GEMINI calculation. Red circles: 
filtered calculated data. The large effect of detector filter on charge Z=4 is due to the unbound 
8 Be isotope that is present in the calculation but filtered because not identified in data analysis. 



cos(9prox) for the primary dynamically emitted fragments (5<Z<8) with the asy-stiff (open 
blue circle) and asy-soft (open red squares) parametrizations, respectively. Most of the 
dynamically emitted fragments are produced at values of cos(Qprox) >0.8 with increasing 
value of <N/Z> by increasing the degree of alignment. It is interesting to note that the asy- 
stiff parametrization tends to produce more neutron rich fragments and with a steeper slope 
with respect to the asy-soft one. Results of SMF+GEMINI are shown in Fig. [7K) as hatched 
zones for asy-stiff (dark-blue) and asy-soft (grey-red) respectively. Notice that for values 
of cos(6'p/ ?O x)>0.9 the signal of neutron enrichment persists also after the de-excitation 
stage. Comparison between data and calculation shows, as in Fig. HI that the asy-stiff 
parametrization matches better the experimental data. These results are consistent with 



recents measurements obtained from isospin diffusion studies 



1CJ, heavy residue production 



in semi-central collisions |49|] and transverse collective flow of light charged particles [12 ] . 
Although final results of <N/Z> observable depend also upon the amount of excitation 
energy given as input in GEMINI calculations, and, consequently, secondary decay tends 



generally to reduce the sensitivity to symmetry energy 
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481 ] . we have shown that it is possible 




FIG. 7. (Color online) a) <N/Z> as a function of cos(9prox) for charges 5<Z<8: solid circles: 
experimental points; SMF calculations are shown for primary fragments as empty circles (asy- 
stiff) and empty squares (asy-soft); SMF-GEMINI calculation for asy-stiff (blue hatched area) and 
asy-soft (red hatched area) parametrization. b) experimental correlation for statistically emitted 
particles (see text). The dashed line gives the average N/Z over the whole cos(9prox) interval. 



to construct robust observable for neck emission dynamics that maintains memory of the 
early stages of the reaction. This work adds new important observable for isospin dynamics 
studies in heavy ion collisions and improves the consistency of the different analyses that 
have been performed so far on the symmetry term of the EOS 
perspectives for reaction studies with exotic beams. 



501 ] . It also opens new 



IV. CONCLUSION 

In summary, we have presented new experimental results correlating the emission 
timescale of intermediate mass fragments (IMF) at mid-rapidity in semi-peripheral col- 
lisions with their isotopic composition. We have shown that large values of <N/Z> are 
acquired by light IMFs dynamically emitted in the early stage of the reaction, for both 



14 



the neutron rich and neutron poor systems studied here. By comparing the experimental 
data with SMF calculations, we have produced valuable information on the parametrization 
of the symmetry energy term of EOS, getting indication for a moderately stiff symmetry 
potential, and new constraints for further simulations of the reaction dynamics. 

ACKNOWLEDGMENTS 

We thank the INFN-LNS accelerator staff for providing beam of excellent quality. We 
acknowledge also C. Marchetta and E. Costa for providing high-quality targets. This work 
for V. Baran was partially supported by the Romanian National Authority for Scientific 
Research, CNCS-UEFISCDI, Project No. PN-II-ID-PCE-201 1-3-0972. 



[1] S. Hudan et al., Phys. Rev. C 71, 054604 (2005). 

[2] S. Piantelli et al, Phys. Rev. C 74, 034609 (2006). 

[3] M. Di Toro, A. Olmi and R. Roy, Eur. Phys. J. A 30, 65 (2006). 

[4] V. Baran, M. Colonna and M. Di Toro, Nucl. Phys. A730, 329 (2004). 

[5] J. Lukasik et al., Phys. Lett. B 566, 76 (2003). 

[6] C.J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 5647 (2001). 
[7] A.S. Botvina and I.N. Mishustin, Nucl. Phys. A843, 98 (2010). 

[8] M.B. Tsang, Y. Zhang, P. Danielewicz, M. Famiano, Z. Li, W.G. Lynch, and A.W. Steiner, 

Phys. Rev. Lett. 102, 122701 (2009) and references therein. 
[9] Z.Y. Sun et al, Phys. Rev. C 82, 051603(R) (2010). 
[10] E. Galichet et al., Phys. Rev. C 79, 064614 (2009); E. Galichet, M. Colonna, B. Borderie and 

M.F. Rivet, Phys. Rev. C. 79, 064615 (2009); 
[11] M.A. Famiano, T. Liu, W.G. Lynch, M. Mocko, A.M. Rogers, M.B. Tsang, M.S. Wallace, 
R.J. Charity, S. Komarov, D.G. Sarantites, L.G. Sobotka, and G. Verde, Phys. Rev. Lett. 97, 
052701 (2006). 
[12] Z. Kohley et al., Phys. Rev. C 83, 044601 (2011). 

[13] J. Rizzo, M. Colonna, M. Di Toro, H.H. Wolter and M. Zielinska-Pfabe, Nucl. Phys. A806, 
79 (2008). 

15 



[14] M. Di Toro, V. Baran, M. Colonna and. V. Greco, J. Phys. G 37, 083101 (2010). 

[15] V. Baran, M. Colonna, M. Di Toro, M. Zielinska-Pfabe, and H.H. Wolter, Phys. Rev. C 72, 

064620 (2005). 
[16] J.F. Dempsey et al., Phys. Rev. C 54 1710 (1996). 

[17] L.G. Sobotka, R.J. Charity, D.K. Agnihotri, W. Gawlikowicz, T.X. Liu, W. Lynch, U. 

Schroder, J. Toke, and H.S. Xu, Phys. Rev. C 62, 031603(R) (2000). 
[18] I. Lombardo et al., Phys. Rev. C 82, 014608 (2010). 
[19] J. Wilczyhski et al., Int. J. Mod. Phys. E14, 353 (2005). 
[20] E. De Filippo et al., Phys. Rev. C 71, 044602 (2005). 

[21] E. De Filippo et al., Proceedings of the IV Workshop on particle correlations and femtoscopy, 

Krakow 2008, Acta Phys. Pol. B 40, 1199 (2009). 
[22] A. Pagano et al., Nucl. Phys A 734, 504 (2004); A. Pagano, Nucl. Phys. News, vol. 22, No 1, 

28-33 (2012). 

[23] N. Le Neindre et al., NIM A 490, 251 (2002); M. Alderighi et al., NIM A 489, 257 (2002). 
[24] C. Cavata, M. Demoulins, J. Gosset, M.C. Lemaire, D. L'Hote, J. Poitou, and O. Valette, 

Phys. Rev. C 42, 1760 (1990). 
[25] P. Russotto et al., Phys. Rev. C 81, 064605 (2010). 

[26] P. Russotto et al., INPC2007 Proceedings (Tokio, June 3-8 2007), edited by S. Nagamiya, T. 

Motobayashi, M. Oka, R.S. Hayano and T. Nagae, pag. 507 (2008). 
[27] D.J. Hinde, J.R. Leigh, J.J.M. Bokhorst, J.O. Newton, R.L. Walsh, J.W. Boldeman, Nucl. 

Phys. A472, 318 (1987). 
[28] M. Papa et al., Phys. Rev. C75, 054616 (2007). 

[29] S. Piantelli, L. Bidini, G. Poggi, M. Bini, G. Casini, P. R. Maurenzig, A. Olmi, G. Pasquali, 

A. A. Stefanini, and N. Taccetti, Phys. Rev. Lett. 88, 052701 (2002). 
[30] S. Hudan, R.T. de Souza and A. Ono, Phys. Rev. C73, 054602 (2006). 
[31] F. Bocage et al., Nucl. Phys. A 676, 391 (2000). 
[32] A.B. Mcintosh et al., Phys. Rev. C 81, 034603 (2010). 
[33] A. Mangiarotti et al., Phys. Rev. Lett. 93, 232701 (2004). 
[34] I. Skwira-Chalot et al., Phys. Rev. Lett. 101, 262701 (2008). 
[35] P.M. Milazzo et al., Nucl. Phys. A703, 466 (2002). 
[36] D. Theriault et al., Phys. Rev. C 71, 014610 (2005). 

16 



[37] R. Planeta et al., Phys. Rev. C 77, 014610 (2008). 

[38] R. Lionti, V. Baran, M. Colonna and M. Di Toro, Phys. Lett. B 625, 33 (2005). 
[39] I. Lombardo et al., Phys. Rev. C 84, 024613 (2011). 

[40] M.V. Ricciardi, A.V. Ignatyuk, A. Kelic, P. Napolitani, F. Rejmund, K.H. Schmidt, O. Yor- 

danov, Nucl. Phys. A 733, 299 (2004). 
[41] E. Geraci et al., Nucl. Phys. A 732, 173 (2004). 
[42] G. Casini et al.. larXiv:1202.Q192i yl (2012). 

[43] E.M. Winchester, J. A. Winger, R. Laforest, E.Martin, E. Ramakrishnan, D.J. Rowland, A. 

Ruangma, and S.J. Yennello, Phys. Rev. C 63, 014601 (2000). 
[44] M. D'Agostino et al., Nucl. Phys. A861, 47 (2011). 

[45] M. Colonna, M. di Toro, A. Guarnera, S. Maccarone, M. Zielinska-Pfabe and H.H. Wolter, 

Nucl. Phys. A642, 449 (1998). 
[46] M. Colonna (private communication). 

[47] R.J. Charity, M.A. McMahan, G.J. Wozniak, R.J. McDonald, L.G. Moretto, D.G. Sarantites, 
L.G. Sobotka, G. Guarino, A. Pantaleo, L. Fiore, A. Gobbi, K.D. Hildenbrand, Nucl. Phys. 
A483, 371 (1988). 

[48] P. Napolitani, M. Colonna, F. Gulminelli, E. Galichet, S. Piantelli, G. Verde, and E. Vient, 

Phys. Rev. C 81, 044619 (2010). 
[49] F. Amorini et al., Phys. Rev. Lett. 102, 112701 (2009); G. Cardella et al., accepted for 

publication in Phys. Rev. C (2012). 
[50] M.B. Tsang et al.. larXiv: 1204.04661 (^2012). 



17 



